Perturbative Effective Theory in an Oscillator Basis? 



W. C. Haxton and T. Luu 
Institute for Nuclear Theory, Box 351550, and Department of Physics, Box 351560, 
University of Washington, Seattle, WA 98195 
(February 8, 2008) 

I 

< 

. The effective interaction/operator problem in nuclear physics is believed to 
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be highly nonperturbative, requiring extended high-momentum spaces for ac- 

t-H ■ 

. curate solution. We trace this to difficulties that arise at both short and long 

CN 
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. distances when the included space is defined in terms of a basis of harmonic 

oscillator Slater determinants. We show, in the simplest case of the deuteron, 
that both difficulties can be circumvented, yielding highly perturbative results 
in the potential even for modest (~ 6Hlo) included spaces. 
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There is an extensive literature on attempts to relate the effective interaction, needed 



in any description of nuclei based on a finite set of low-momentum basis states, to the 
underlying, more singular iViV interaction. As Barrett |1[ has summarized, one hope of 
investigators in the 1970s was that H e ^ might be expanded perturbatively in either the bare 
potential V or in the G matrix, the sum of all two-nucleon ladder diagrams for scattering in 
the excluded, high-momentum space. While some phenomenological success was achieved 
by selecting certain diagrams @, more systematic treatments provided little indication of 
convergence. For example, the third-order calculations of Barrett and Kirson || yielded 
a correction of the same size but opposite in sign to the second-order result. At about 
the same time Shucan and Weidenmuller [Q demonstrated in a toy model that strongly 
coupled "intruder" states - states low in energy having little overlap with the model space 
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- generically lead to poorly convergent expansions. These disappointing results led many 
practioners to turn to phenomenological H e ^s - an approach that is also lacking because it 
fails to provide a basis for generating effective operators consistent with H e ^ . 

The lack of a first-principles technique for effective interactions had a discouraging effect 
on the field for a number of years. However several recent developments - including the 
rapid growth of computing power and interest in effective field theory || - have encouraged 
new attempts to calculate H e ^ in a systematic, controlled way. 

One step in this direction was reported recently ||: for an included space consisting of 
a complete set of harmonic oscillator (HO) Slater determinants of energy < Asm^, the 
deuteron and 3 He were treated as exact effective theories (ET). The point was to illustrate 
crucial aspects of ETs that are generally absent in models like the shell model (SM), including 
nontrivial wave function normalizations, the many-body nature of effective interactions, the 
rapid evolution of H e ^ matrix elements with changes in the model space, and the crucial 
role of consistent effective operators. 

The starting point for the calculation of Ref. is the Bloch- Horowitz equation JF]. For 
a low- momentum "included space" defined by Asm and a HO size parameter b, H e ^ will 
be translationally invariant. (Note %u = h 2 /Mb 2 , with M the nucleon mass.) Defining the 
projection operator onto the high-momentum Slater determinants by Q(AsM,b), one finds 

H "' = H + H E^QH QH 

H'"\4>sm) = E\^ SM ) \4, SM ) = (1 - Q)W) (1) 

where is the exact wave function, H\ip) = E\ip), and H = T + V is the sum of rela- 
tive kinetic energy and potential terms. This equation has to be solved self-consistently for 
each desired but unknown eigenvalue E. The method introduced in |J provided an efficient 
solution to the self-consistency problem, constructing the needed Green's function as a func- 
tion of E by a method based on the Lanczos algorithm: with this technique one iteratively 
extracts the high-momentum spectral information most relevant to the construction of the 
Green's function. 
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This approach, in principle, can be extended to heavier systems. In practice, however, 
the integration over high-momentum states becomes more challenging with increasing A: in 
||, where the Argonne v 18 potential || was employed, high- momentum integrations up to 
Aqo ~ 140 (~ 60) were necessary to achieve 1 keV (25 keV) accuracy. In this letter we argue 
that this integration can be performed in a much simpler way. 

We begin by investigating why the high momentum integration is nonperturbative. One 
can envision moving the SM scale Asm very close to the necessary Aoo, which clearly makes 
the excluded-space contribution to the BH equation a small correction. As the energy 
denominator in Eq. (1) is then very large, one might assume that the high-momentum 
problem is now perturbative. 

To test this conclusion we expand the excluded-space contribution to Eq. (1) 



H QH = H 

E-QH^ 



1 + —Jl—QV—!—— + 



QH (2) 



E-QT E-QT^ E-QT 

The order-by-order results ( dashed lines in Fig. la) are quite curious. The total excluded- 
space contribution, given the very high value of Asm — 70 chosen, accounts for only 20 keV 
of binding energy; 85% of this is generated by the first term in Eq. (2). However subsequent 
order-by-order corrections appear to converge to a value a few keV above the true binding 
energy. Only after ~ 1000 orders of perturbation theory is the correct binding energy slowly 
achieved. (Very similar results are obtained if one uses the HO Hamiltonian T + Vq, instead 
of T, as the unperturbed Hamiltonian.) 

In Figs, lb and lc this behavior is explored for two matrix elements of H e ^ . Most 
matrix elements converge rapidly, like the Is — Is example of Fig. lb. The exceptions are 
those where the bra or ket lies in the last included shell, i.e., A a = Asm, such as the n = 36 
s-state case of Fig lc). This associates the poor convergence with T which, because of the 
raising/lowering properties of the gradient operator in a HO basis, connects only such states 
to the excluded space. 

It is helpful to note that a HO-basis ET differs from effective field theories in that the 
expansion is around an intermediate momentum q ~ 1/6, rather than q ~ 0. Fig. 1 



shows that this expansion then induces a familiar long-distance problem: because HO wave 
functions fall off too sharply, the correct asymptotic wave function can only be achieved by 
scattering through a very large number of high-momentum oscillator states. Consequently a 
poorly convergent "tail" exists regardless of how high Asm is chosen: a larger Asm restricts 
the unresolved tail to larger r and thus limits its numerical significance, but does not in 
principle make the problem perturbative. 

The solution to this problem is not entirely trivial. The HO basis is essential because 
of the center-of-mass separability it provides. Because the relevant operator is the relative 
kinetic energy 

T=\jlT l3 (3) 

where T tj = (p* - Pj) 2 /2AM, the missing long-distance correlations are two-body. As the 
problem is associated with T hopping from the included space at large r with a large am- 
plitude, we rearrange the BH equation so that the Green's function is sandwiched between 
Vs, thereby cutting off large-r propagation || 

(a\H*"\fl = (a\T\P) + (a\V + V R -L^ QV\fi (4) 

where 

^ ) = E—Qf^ (5) 
with \a) a HO Slater determinant. Thus |a) = \a) for Slater determinants for which A a ^ 
Asm- Otherwise, apart from the overall normalization, \a) differs significantly from \a) only 
in its large-r behavior (Fig. 2). We stress that Eq. (4) is equivalent to Eq. (1). Finally we 
insert the bracketed Green's function expansion of Eq. (2) into Eq. (4). Thus QV always 
appears between insertions of QT, summed to all orders. 

The results are given by the solid lines in Figs, la and lc. The (36s\H e ^\36s) matrix 
element and total binding energy now converge rapidly. In fact Asm can be lowered to ~ 40 
while keeping high momentum contributions perturbative: in third order 1 keV accuracy is 
maintained. 



However, with further lowering of the SM scale, the convergence again exhibits hints of 
deterioration - with new symptoms. For A S m = 30, the third-order calculation reduces ~ 
10% errors in the bare value of (ls\H e ^\ls) to ~ 0.2%. But an error in excess of 0.1% - 
corresponding to 50 keV - persists after 10 additional orders of perturbation. Unlike the case 
discussed above, all matrix elements are affected, though the nonperturbative corrections to 
s — s transitions are larger than those for s — d and much larger than those for d — d matrix 
elements. This is a signature of scattering at small r through the hard core. Numerically 
one can verify that the nonperturbative tail disappears if the hard core in V is removed. 

An exact ET must yield the same result regardless of the choice of excluded space 
parameters b and A SM . However it is possible that a judicious choice of these parameters 
might simplify the numerical difficulty of the ET. While A S m generally is limited by the size 
of one's computer, b can be varied freely. A natural choice for b is the value that minimizes 
the 0th order energy (obtained by ignoring entirely the last term in Eq. (1) or Eq. (4)). This 
corresponds to minimizing the bare Hamiltonian (T + V)sm as a function of b in Eq. (1). 
The closer the binding energy to the correct value, of course, the smaller the contribution 
of the high-momentum corrections due to scattering in the excluded space. 

The dashed lines in Fig. 3 are the Oth-order energies for Eq. (1) as a function of b 
for Asm = 6, 8, and 10. As (T + V)sm is a variational estimate of the energy, the contour 
minima are upper bounds to the exact result. The Oth-order results fail to bind the deuteron; 
the minima are achieved for b ~ 0.79-0.83/. (Note that b e ff = \^2b is the size scale in the 
relative coordinate, f[ — f^.) This fe e // is considerably smaller than the nuclear size: the 
SM is doing its best to find a compromise between two needs, resolving the hard core (a 
problem that becomes easier for small b) and reproducing the correct long-distance behavior 
(a problem that becomes easier for large b - a doubling of b roughly halves the number of 
high-momentum states that must be included to calculate T e ^ to an equivalent accuracy). 
The resulting compromise b addresses neither need well. 

The solid lines in Fig. 3 are the corresponding results for the Oth-order approximation to 
Eq. (4). The minima shift to b ~ 0.4-0.5/, and the unperturbed results are very accurate, 
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with errors of 21, 36, and 52 keV for 10, 8, and 6 fuv spaces - an improvement of about a 
factor of 100 in binding-energy accuracy over the corresponding dashed-line results. This 
has been achieved without taking into account any effects of QV. The interpretation is 
clear: once one has solved the long-distance problem through the resummation of Eq. (4), 
the nonperturbative effects (as we shall see below) of the hard core can be absorbed almost 
entirely into the included space, given an appropriate HO size scale b ~ r c , where r c is the 
hard-core radius. This scale arises naturally out of the minimization, because the short-range 
repulsion is such a dominant feature of the potential. 

The Ahu Oth-order energy minimum is not quite so impressive: the included space is 
sufficiently restrictive that a nonnegligible contribution from QV remains even if b is opti- 
mized. Conversely, the minimum at lOhu is very flat, which means that there is a range of 
b - a set of included spaces - in which the effects of QV remain very small. 

Fig. 4 gives our most important results. Fig. 4a illustrates the problem noted in the 
introduction, a perturbative expansion of Eq. (1) converges very poorly. Only at 10 Tiuj 
(not shown) do the fluctuations begin to damp out, though they remain at the ~ 100 keV 
level even at high orders. In contrast, the expansion based on Eq. (4) (Fig. 4b) is highly 
perturbative. The lst-order correction V(E — QT)~ X QV for the 6,8, and 10 fwj calculations 
yields binding energies accurate to <3 keV. We see that the 10 and 8 Tiuj calculations are 
effectively identical at and beyond 2nd order; the 10, 8, and 6 Tiuj results merge at 4th order. 
Even the ATtuj Oth-order result is quickly corrected to an accuracy of 1.2 keV at 3rd order. 

The accuracy of the unperturbed lOhu result of Eq. (4) is comparable to that achieved 
in || by direct summation of the high-momentum contribution to Eq. (1) to A ro ~ 70. 
Perhaps more important, however, is the promising result that the remaining corrections 
associated with QV are perturbative. This has important implications for the complexity 
of H e ff in heavier nuclei, suggesting that the number of nucleons in the excluded space 
that must be linked by QV can be limited. For example, V(E — QT)~ X QV links only three 
nucleons. Rapid convergence in the perturbation expansion translates into a hierarchy of 
three-body, four-body, etc, contributions to V e ^ of rapidly diminishing importance. 
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So far the discussion has emphasized results and their interpretation rather than technical 
aspects. However, the technical aspects deserve discussion because many HO properties can 
be usefully exploited: 

1) The operator T is diagonal in a plane wave basis. The transformation into that basis is 
particularly simple for the HO, as the Fourier transform of a HO is a HO. 

2) The operator appearing in the Green's function, however, is QT, which is more difficult 
to treat. One can relate this Green's function to the free Green's function G = ^z^, at 
the cost of a matrix inversion in the included space. Using the included-space projector 
P = 1 — Q we define operators acting within the included space 

T = PG P, r\ = PG VG P, T 2 = PGoVGoVGoP ... (6) 

One can then write the perturbative expansion of H e ^ as 

ir" — (e — r l ) + r^ro 1 + rv 1 ^ - r^- 1 

+ rv 1 ^ - r^o 1 ^ - r 2 r 1 r 1 + r^- ^rv 1 + ... (7) 

where the terms correspond to the contributions from T, V, (QV) 2 , (QV) 3 , respectively. 
Note that the matrix Tq 1 differs from the matrix P(E — T)P only in the entries where the 
bra and ket belong to the last included "shell" . 

3) As T is a sum over all relative momenta, an evaluation of G in an independent-particle 
basis will generate a product of two-particle correlations (with a dependence on each relative 
partial wave). In contrast, with Jacobi coordinates G is diagonal in momentum space 

G ° = E-^{kl + ... + k\_ l ) (8) 

where hi is the momentum associated with the ith Jacobi coordinate. Clearly the Jacobi 
basis is the simpler choice. As noted before, if we work to some order in perturbation 
theory, then the number of nucleons interacting at one time in the excluded space will 
be limited. Thus, in general, there will be "spectator" nucleons, and the matrix element 
will have a spectator dependence corresponding to the overlap integral weighted by r . The 
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importance of such spectator dependence in a HO-based ET was noted in |1(J • The integrals 
over spectator nucleons can be handled analytically and are reducible to finite sums over 
gamma and incomplete gamma functions. 

These properties allow one to evaluate the HO matrix elements of the operators of Eq. 
(7), and thus the binding energy term by term in perturbation theory. In the present case of 



the deuteron the relative-coordinate matrix element is transformed into momentum space 



(n f (£fs)jmtm t \T a (E)\rii (£is)jmtm t ) 



hu eff 



4(n/-l)!(ni-l)! 



r(n f +£ f + l/2)T(n i +£ i + l/2) 



1/2 



3C 



x / dp f Pf 6 L e n f+ }f(p 2 f ) [ dpi Pl 6 L l ^}l 2 (p^){p f (£ f s)imtm t \V a \pi(£is)imtm t ) (9) 
Jo pj + E J Jo px + E 

where u e ff = uj/2, L is a Laguerre polynomial, and Eq is the dimensionless binding energy 
\E\/hu e ff. The momentum-space matrix element of the iterated potential is given by a 
simple recursion relation 

roc p 2 

{p f {£fs)jmtm t \V a \p i (£ i s)jmtm t ) - V / dp {p f {£ f s)jmtm t \V 1 \p{£s)jmtm t ) 

- JO P Z + tlQ 

x (p(£s)jmtm t \V a „i\pi(£is)jmtm t ) (10) 

where (p f(£ fs)jmtmt\Vi\p(£s)jmtmt) is given by 

1 f°° I~2 l~2 

/ drr 2 \ -j £ (p f r)((£ f s)jmtm t \V(rb eff )\(£s)jmtm t )\ -jg(pr) (11) 

nuJeff Jo V 71 V 71 

Note that all integration variables are dimensionless. 

In summary we have shown that conventional HO-basis effective interactions calculations 
involve both long- and short-distance nonperturbative scattering. The effects of such scat- 
tering can be absorbed into the included space by an appropriate summation of the relative 
kinetic energy to all orders, followed by a tuning of the HO included space to absorb most 
of the hard-core scattering. This tuning results from minimizing the Oth-order energy as a 
function of b. In the test case of the deuteron, we find very accurate 0th order results and 
very rapid convergence in further orders of perturbation. Existing Jacobi coordinate SM 



codes can treat light lp-shell nuclei and three-body interactions [11]], which should allow 



similar calculations to be done through lst-order in the two-body potential (Oth order in 
the weaker three-body potential) in these cases. 

We thank Bob Wiringa for providing the v 18 potential codes, and Silas Beane and Martin 
Savage for helpful discussions. This work was supported in part by US Department of Energy 
grants DE-FG03-00ER41132 and DE-FC02-01ER41187. 
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FIG. 1. The dash lines were obtained by inserting the Green's function expansion of Eq. (2) 
into Eq. (1). Very slow convergence is found for the binding energy (a) and for the effective matrix 
elements that involve included states in the last "shell." Other matrix elements (b) converge 
quickly. The corresponding results for Eq. (4) (solid lines) all converge rapidly. The horizontal 
lines are the exact results. 



11 




1 2 3 4 5 6 7 



r(fm) 



FIG. 2. For A 5M = 10 the last included s-state has n = 6. The state E/(E - QT)\n = Q£ = 0) 
(solid line) is compared to its HO counterpart, \n = 6£ = 0) (dashed line), showing that the former 
is modified at large distances. The modified wave function was multiplied by 1.934 to match the 
HO at r = 0. 
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FIG. 3. The 0th order approximations to Eq. (1) (dashed lines) and Eq. (4) (solid lines) as a 
function of b, showing the shifts in the minima and improved binding energies associated with the 
latter. 
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FIG. 4. Binding energies obtained by treating the excluded-space contributions to Eq. (1) (a) 
and Eq. (4) (b) perturbatively, via the Green's function expansion of Eq. (2). The solid, dotted, 
dashed, and dot-dashed lines correspond to 10, 8, 6, and 4 hco included spaces, respectively. 
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